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ABSTRACT 

Deep infrared and optical images are presented of three proposed remnants of 
Thorne-Zytkow Objects (TZO). In particular, the IR data go several infrared magni- 
tudes deeper than previous observations and in at least one case reveal the existence 
of weak objects within the error circles. It is argued, however, that none of the objects 
is likely to be the binary companion to the X-ray source in that region. These data 
present severe limits on any possible star or residual envelope at the distances of the 
respective X-ray objects. 

Key words: circumstellar matter - infrared: stars. 



1 INTRODUCTION 

Recent work by several authors, especially van Paradijs et 
al. (1995) and Mereghetti & Stella (1995), have identified a 
group of pulsating X-ray sources that appear to have no de- 
tectable optical counterpart. In addition, these objects show 
no evidence of any binarity and all have pulse periods in the 
5- 10s range. Their X-ray luminosities are all sufficiently high 
that accretion must be occurring on to the neutron star from 
some source other than the interstellar medium. One possi- 
ble scenario is that they are the first objects identifiable as 
the remains of Thorne-Zytkow Objects (TZOs). 

The progenitor TZOs would have consisted of a neutron 
star embedded in the core of another star (Thorne &Zytkow 
1977). It is thought that these TZOs could be formed as 
a result of the evolution of high mass systems through a 
common envelope phase (Ghosh et al. 1997). Once formed, 
thermonuclear burning via the r-p process occurs in the core 
region of the TZO. This supports the highly convective outer 
envelope against gravitational forces. However, once there is 
insufficient fuel remaining for the r-p process to occur effi- 
ciently, the outer envelope collapses on the Kelvin-Helmholtz 
timescale (Podsiadlowski, Cannon and Rees, 1995). Some of 
the envelope material then forms into a disk around the neu- 
tron star with a radius that is dependent on the amount of 
angular momentum in the system. It is from this disk of 
material that the neutron star may accrete enough matter 
to power the X-ray pulsations seen. 

In order to explore this possibility deep infrared imag- 
ing observations have been carried out from UKIRT and the 
NASA IRTF to search for any IR emission from such a disk 
around three of the neutron stars. The observations reach as 
deep as J~20 and in some cases reveal the presence of faint 



objects in the X-ray error circles. However, the character- 
istics of the IR colours, and other considerations, strongly 
suggest that these are not the TZO envelopes around the 
neutron stars. However, the new stringent limits on any IR 
companion in these systems leave these objects as enigmatic 
as ever. 

The three systems investigated here are: 



1.1 1E2259+586 

The X-ray pulsar 1E2259+586 lies in the centre of the su- 
pernova remnat G109. 1-1.0. The most recent estimates on 
its age make it approximately 3000 years old (Parmar et al. 
1997), though values up to 10000 years have been quoted. 
G109. 1-1.0 is at an estimated distance of about 4 kpc. 
Fahlmanet al. (1982) demonstrated that the X-ray column 
densities to both the pulsar and the SNR were similar and 
hence that these two objects were probably related. How- 
ever, the latest results from Parmar et al.(1997) indicate 
that this may not be the case. They find the X-ray column 
density to the pulsar to be greater than that to the super- 
nova remnant. They suggest that this could be due to the 
pulsar being at a greater distance than the supernova rem- 
nant or to the presence of absorbing material local to the 
pulsar. 

However, assuming the same distance to 1E2259+586 as 
the SNR implies a pulsar X-ray luminosity of w 10 28 Watts. 
Results from EXOSAT, GINGA and BeppoSAX (Hanson et 
al. 1988, Koyama et al. 1989 and Parmar et al. 1997 re- 
spectively) have shown that the 7s period of 1E2259+586 is 
slowly increasing, but the spindown is too slow to power the 
observed luminosity. There is no evidence of binary motion, 
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the upper limit on a x sin i from pulse timing measurements 
being as low as 30 light-ms (Mereghetti, Israel and Stella, 
1998). Early searches for the optical counterpart (Fahlman 
et al. 1982; Margon & Anderson 1983) found several possible 
candidates, one of which (star D) was tentatively identified 
with 1E2259+586 by Middleditch, Pennypacker & Burns 
(1983) on the basis of IR pulsations. However, more sensitive 
fast photometric observations by Davies et al. (1989) failed 
to find any pulsations, and there would now appear to be 
little positive evidence associating star D with 1E2259+586. 
A fainter, multicolour photometric study of all the possible 
counterparts by Davies & Coe (1991), together with revised 
astrometry, updated the list of possible counterparts. 

Coe & Jones (1992) presented the first optical spec- 
troscopy of all the candidates brighter than V=23 and re- 
analysed the Einstein X-ray images of this source taken in 
1981 in order to check the position and size of the X-ray error 
circle. They found that the relocation and reduction in the 
size of the error circle stemming from the re-processing of 
these data considerably changes the perspective on possible 
counterparts to the X-ray pulsar. In addition, a careful study 
of the images reveals previously unreported features that are 
possibly associated with jet activity from the pulsar, a possi- 
bility originally suggested by Gregory and Fahlman (1980). 
However, Hurford & Fesen 1995 report subsequent Rosat 
HRI observations which, while confirming the existence of 
such features, they suggest are just statistical fluctuations 
in the SNR emission. 



1.2 4U0142+62 

This persistent source has appeared in most X-ray cata- 
logues and was originally confused (due to poor angular 
resolution) with the Rosat source RX J0146. 9+6121 (=LSI 
+61 235) - the two sources are only 24 arcminutes apart. 
Measurements by White et al. (1987) with EXOSAT re- 
ported a 25m modulation from the region (later shown to 
be coming from LSI +61 235 (Mereghetti, Stella & De Nile 
1993)) but using the unmodulated lower energy data they 
were able to identify the position of 4U0142+62. The de- 
tection of X-ray pulsations at 8.7s was reported by Israel et 
al. (1994) who concluded that they were most likely com- 
ing from the optically unidentified system 4U0142+62 rather 
than the new Rosat source. Further measurements by White 
et al. (1996) using ASCA confirmed this conclusion and 
also demonstrated evidence for an excess in the X-ray halo 
around the source. This could be due to material local to 
the pulsar along the line of site, possibly in the form of a 
molecular cloud or perhaps the remains of a TZO envelope. 
No evidence of any binary motion has been reported from 
the X-ray signal. 

The accurate position reported by White et al. (1987) 
allowed them to locate the X-ray error circles from both EX- 
OSAT and Einstein on the sky. As in the case of 1E2259+586 
there is no evidence of any obvious optical candidates within 
the defined regions. They obtained an R band image and 
set a limit of R>22.5 for any counterpart. The ASCA ob- 
servations (White et al. 1996) determined the column to the 
source to be 8 x 10 21 cm" 2 which in turn allows them to 
deduce an optical extinction of Ai,~4.7. The distance to the 
source is not well defined but probably lies in the range 0.5- 
2.0 kpc. As in the case of 1E2259+586 a molecular cloud 



N 




Figure 1. A K band image of the field of 1E2259+586. The field 
size is 73 arcseconds and the circle shown represents the Rosat 
X-ray error circle of radius 5 arcsec (Dennerl, private communi- 
cation). The three very faint objects within this error circle are 
referred to in the text as objects L, K and N from East to West 
respectively. 



lies near to, or in front of 4U0142+62 (Mereghetti & Stella 
1995) and may be affecting the column measurement. 



1.3 RX J1838. 4-0301 

This source was originally reported by Schwentker (1994) as 
5.5s pulsar possibly associated with a SNR at a distance of 
4kpc. Interestingly, in this case the X-ray positional error 
circle was reported to include a ~14 mag optical candidate. 
Subsequently Mereghetti et al. (1997) reanalysed the same 
Rosat data and reported are different interpretation of the 
data. They obtained an optical spectrum of the star in the 
error circle and concluded that this object was simply a K5 
main sequence star and that the X-ray emission originated 
from its corona . 



2 THE OBSERVATIONS 
2.1 1E2259+586 

This source was observed by UKIRT using IRCAM3 on 1996 
October 23 in Service mode. The plate scale was ~0.3 arcsec 
per pixel giving an image size of 73 x 73 arcsec. Observa- 
tions were taken in both J and K bands using a standard 
mosaicing technique with a total exposure of 27 minutes in 
each band. The K band image of the field is shown in Figure 
1. 

Superimposed on the image is the Rosat error circle 
(Dennerl, private communication) and the labels for several 
of the brighter stars in the region (following the convention 
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Table 1 . Infrared observational measurements of the candidates 
for 1E2259+586 



Object 


J band 


K band 


(J-K) 


X 


17.0±0.1 


16.1±0.06 


0.89±0.12 


w 


18.1±0.1 


17.0±0.1 


1.1±0.14 


L 


20.1±0.5 


17.9±0.4 


2.2±0.6 


K 


20.0±0.5 


18.0±0.4 


2.0±0.6 


N 


>20.6 


18.3±0.4 


>2.3 



of Coe et al. 1994a). The image is approximately 1-2 mag- 
nitudes deeper than any previous measurements and for the 
first time permits measurements of the three faint objects 
L, K and N. 

In order to obtain an accurate calibration of the images, 
the established objects Rl, R2, G, F, A and D (all shown in 
the figure) were used. The instrument magnitudes of these 
objects were compared with the reported fluxes in Davies & 
Coe (1991) and an excellent correlation was obtained (this, 
incidentally, confirms that there has been no detectable vari- 
ation in the fluxes from any of these objects). Using this 
correlation the magnitudes for the three candidates L, K 
and N, as well as the nearby objects X and W, were deter- 
mined. The results are presented in Table 1. From sampling 
the background around the area of interest the three sigma 
upper limits for any undetected source in this field may be 
determined to be J>19.6 and K>18.4. 

2.2 4U0142+62 

This source was observed by the NASA IRTF on 1996 Au- 
gust 25 using the NSFCAM imager. This instrument has a 
plate scale of 0.3 arcsec per pixel giving an image size of 
77 arcsec. Observations were taken in the J and K bands 
of 30 minutes each. In addition to the target field the stan- 
dard star HD18881 was observed together with an appropi- 
ate number of flats. 

The K band image of the field is shown in Figure 2. Also 
shown in this figure is the Rosat positional error circle from 
Hellier (1994). As one can see the field is extremely empty 
lacking any possible candidates close to the X-ray position. 
The faintest object detected and measured had a magnitude 
of J=19.6±0.1 and K=16.88±0.07. The measurements of all 
the objects identified in the figure are given in Table 2 (note 
that star identification scheme is an extension of that started 
by White et al 1987). 

2.3 RX J1838. 4-0301 

Notwithstanding the discussion in Section 1 of this paper, 
this source lacks a securely identified counterpart. Even the 
question of the existence, or otherwise, of X-ray pulsations 
needs resolving. It was therefore felt appropiate to explore 
the region around the error circle to see if any sources re- 
vealed themself as typical of a counterpart to an X-ray 
source. Consequently, on 1996 November 15 images were 
taken in each of the B,V,R and I band using the Tek8 CCD 
camera on the lm SAAO telescope with standard Johnson 
filters. In addition an image was recorded through the Ha 
filter. The exposure times for B,V and I were 300s, while R 
was 200s and Ha was 900s. 
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Figure 2. The upper panel shows a K band image of the field 
of 4U0142+62. The field size illustrated is 71 arcseconds and the 
circle shown represents the Rosat X-ray error circle of radius 3.2 
arcsec (White et al. 1987). The lower panel is a V band image 
(Stcinle, private communication) indicating the identification of 
stars in the region of the X-ray error circle. 

If the X-ray pulsations are real then one obvious expla- 
nation for this system is that it is either a Be/X-ray binary 
system or a supergiant system. If this is the case then the 
optical counterpart would be a strong source of Ha emision. 
To investigate this possibility all the objects on the frame 
were measured in both the R band and Ha and the two 
sets of instrument magnitudes plotted against each other. 
This approach has been very successful in identifying coun- 
terparts to high mass X-ray binaries (see, for example, Coe 
et al., 1994b). In this case the result revealed no obvious 
source of Ha emission in the field. 
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Table 2. Infrared observational measurements of the candidates 
for 4U0142+62 



Table 3. Infrared observational measurements of the candidates 
for RX J1838.4-0301 



Object 




K bciiid 




1 


1 5m +0 005 


1 3 1 26+0 0(14 


375+0 006 


2 


1 5 33+0 01 


14 01 +0 01 


418±0 015 


2 a 


17(11 +0 04 


i 7 nq +0 03 


n 5i ^+o 05 


3 


1 5 on +n 02 


1 5 1 r _|_n ni 


n vfi _|_n f|2 




1 cr 1 /l-Un 01 


i a Qf. _i_n m 


n 97 -Ufl m 


6 


14.95±0.01 


14.63 ±0.01 


0.32 ±0.01 


8 


15.74±0.01 


15.44 ±0.01 


0.30 ±0.02 


9 


15.52±0.01 


15.20 ±0.01 


0.32 ±0.02 


13 


16.82±0.02 


16.38 ±0.02 


0.43 ±0.03 


14 


17.15±0.03 


16.67 ±0.03 


0.47 ±0.04 


16 


17.02±0.03 


16.38 ±0.02 


0.64 ±0.03 


25 


18.55±0.06 


17.54 ±0.04 


1.01 ±0.07 


29 


17.05±0.03 


16.65 ±0.02 


0.40 ±0.04 


30 


17.16±0.03 


16.34 ±0.02 


0.82 ±0.03 




Figure 3. K band image of RX J1838. 4-0301. The field size is 73 
arcseconds and the small circle of radius 5" shows the estimate 
of the location of the Rosat source (Mereghetti et al., private 
communication and 1997), whereas the larger circle (radius 10") 
is that of Schwentker (1994). 



In addition to the above, deep IR observations were ob- 
tained from UKIRT using IRCAM3 on 1997 September 13 
(J band) and 1997 September 29 (K band). Total exposures 
of 18 minutes were used in both bands. The instrument set 
up was identical to the observations of 1E2259±586. The K 
band image is shown in Figure 3 along with the determi- 
nations of the position of the Rosat source by Schwentker 
(1994) and Mereghetti et al. (private communication and 
1997). The J and K magnitudes were measured of all the 
identified objects shown in Figure 3 and the results are pre- 
sented in Table 3. 

It is obvious from the UKIRT image that the reported 
"14th magnitude star" is in fact a complex of 5 or more 
previously unresolved stars. The limits on any undetected 



Object 


J band 


K band 


(J-K) 


A 


12.69±0.03 


12.74±0.03 


-0.05±0.04 


B 


12.29±0.03 


12.62±0.03 


-0.33±0.04 


c 


14 01 +0 03 


14 69±0 03 


-0 68+0 04 


D 


12.03±0.03 


11.98±0.03 


+0.05+0.04 


E 


13.72±0.03 


14.16±0.03 


-0.44±0.04 


F 


13.38±0.03 


13.63±0.03 


-0.25±0.04 


G 


13.36±0.03 


13.24±0.03 


+0.12+0.04 


H 


12.17±0.03 


12.24±0.03 


-0.07+0.04 


.1 


11.55±0.03 


12.01±0.03 


-0.46±0.04 


K 


12.47±0.03 


12.19±0.03 


+0.28+0.04 


L 


13.19±0.03 


14.15±0.03 


-0.96+0.04 


M 


13.61±0.03 


14.24±0.03 


-0.63+0.04 



sources in the field were determined by sampling the back- 
ground around the X-ray error circle. The resulting three 
sigma upper limits are J>17.7 and K>16.9. 



3 DISCUSSION 

In order to investigate the nature of the objects in the 
fields of our targets an infrared colour plot was produced 
for each region. Figure 4 shows the plot for the first source, 
1E2259+586. 

To help interpret this diagram the location of the Main 
Sequence (MS) has been determined with the appropiate 
reddening. The solid line in the figure shows the MS with a 
reddening of A v —2.5 (corresponding to a distance of 4kpc), 
while the dotted line is for A„=4.7 (corresponding to lOkpc). 
These two values for reddening encompass the range of dis- 
tances proposed for 1E2259+586, though the smaller value 
is the more probable (Davies & Coe 1991). As one can see, 
the reddening approximately follows the direction of the MS 
and hence does not provide a very useful distance discrim- 
inator. However, it is clear that most of the more obvious 
sources lie comfortably along the MS line and present no un- 
usual IR colours. Object E, discussed in Coe & Jones (1992) 
as a possible white dwarf, is not detected in either IR image 
and hence is fainter than J=21 and K=18.5. 

The same cannot be said though for the new objects 
labelled K, L and N - these are the objects which actually lie 
in the Rosat error circle. A very large amount of reddening is 
required to make these objects consistent with a MS object - 
for example, a cluster of young B0 stars would need A„~13- 
14. Though it is highly improbable that the reddening in this 
galactic direction could, in general, be larger than A„=~5, 
such a high value could be consistent with a large amount 
of local material around a distant cluster. 

Alternatively, A„=4.7 corresponds to E(J-K)=0.8 
(Rieke & Lebofsky 1985) and the observed colour would 
be consistent with a Seyfert galaxy which typically has 
(J-K) =1.6, and hence would be reddenned to (J-K)=2.4. 
However, the chance of three such objects lying within 10 
arcsec seems unlikely. 

For the field containing 4U0142+61 the brightest 14 
objects have had both their J and K colours determined. 
The results are presented in Figure 5. 

Various fits to the data were undertaken to obtain the 
best estimate of the distance and reddening to the field 
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Figure 4. The infrared colours (J-K) for all the objects close to, 
or in, the X-ray error circle of 1E2259+586. The solid and dotted 
lines represent the location of the Main Sequence at reddenings 
of A„=2.5 and 4.7 respectively. 



stars. The solid line in the figure shows the best estimate 
for the Main Sequence at the distance (lkpc) and reddening 
(A„=4.7) quoted in White et al (1996). Clearly the field stars 
do not lie anywhere close to this position, nor does simply re- 
ducing the reddening to, for example, A v =2 (analagous per- 
haps to a position just in front of a molecular cloud reported 
by Mereghetti & Stella 1995). The best fit was obtained for 
an A„=2.4 at distance of 6kpc - indicated in Figure 5 by the 
dashed line. From this one can conclude that none of the 
visible stars are related to the X-ray source - perhaps not 
too surprising since they all lie well away from the X-ray 
error circle. 

In the case of RX J1838.4-0301 the lack of any obvi- 
ous Ha source in the immediate region of the X-ray posi- 
tion could be interpretted as support for the conclusions of 
Mereghetti et al. (1997) that the X-ray signal is coronal in 
nature. However, the IR images reported here clearly show 
the presence of several comparably bright stars in the region 
of the Rosat source. Figure 6 shows a similar plot to Figures 
4 and 5, but for the objects in the field of RX J1838.4-0301. 

Schwentker (1994) reports a hydrogen column to this 
source in the range (1-3) x 10 21 atoms per sq.cm. If we use 
this value with the relationship given in Ryter, Cesarsky & 
Audouze (1975) we can deduce an E(B-V) in the range 0.2 - 




.1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 1.1 1.2 1.3 



(J-K) 

Figure 5. The infrared colours (J-K) for all the objects close 
to the X-ray error circle of 4U0142+62. The solid line represent 
the location of the Main Sequence at a reddening of A„=4.7 and 
a distance of lkpc, and the dotted line the MS with A V =2A at 
distance of 6kpc. 

0.4 and hence an E(J-K) in the range 0.1 - 0.2. The locus of 
the Main Sequence at an average reddening of E(,I-K)=0.15 
(i.e. E(B-V)=0.27) and the distance of lkpc suggested by 
Schwentker (1994) is marked by the dashed line in Figure 
6. Many of the stars in the region of the X-ray source lie 
close to this line, but the spread is quite significant in the 
(J-K) axes. Mereghetti et al (1997) have proposed a K5 star 
at 300pc as being consistent with the X-ray data (if the 
pulsations are not real). A study of all the stars in this di- 
rection indicates that at distances beyond 200-300 pc the 
hydrogen column will always be of the order 10 21 atoms per 
sq.cm. The solid line in Figure 6 shows the reddening ex- 
pected for such a column, but now at a distance of 300pc. 
Star A lies extremely close to the position expected for a 
K5 star at such a distance and reddening, and this object 
may have been the dominating star in Mereghetti et al's 
data. However, the other stars in the Rosat error circle lie 
scattered all over the diagram and offer many other possi- 
ble solutions. For examply, Star D is the brightest in our IR 
data at J=12.0 and probably lies even closer than Star A. 
Certainly Mereghetti et al might well have had several of 
the stars A-E in the 2" slit for their spectroscopy. 

If the pulsations turn out to be real as proposed by 
Schwentker (1994), then it would support RX J1838.4-0301 
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Figure 6. The infrared colours (J-K) for all the objects close to, 
or in, the X-ray error circle of RX J1838. 4-0301. The solid line 
represent the location of the Main Sequence at a reddening of 
A„=0.42 and a distance of 0.3kpc. The dashed line is the MS at 
the same A^, but at a distance of l.Okpc. 



remaining a member of this class of enigmatic objects. What 
is clear is that it is very unlikely to be a member of any 
Be or supergiant binary system. Further X-ray observations 
are clearly required both to refine the X-ray position in this 
crowded field, and to clarify the situation regarding the pos- 
sible X-ray pulsations. 

Finally, a rough estimate of the limiting size for a disk 
producing black body emission may be obtained from the IR 
results presented here. In the case of 1E2259+586 the limits 
of J>20 and K>18.5 imply a limiting size of approximately 
1 solar radius (7 x 10 8 m) for a 2500K black body at 4.5kpc. 
Podsiadlowski, Cannon and Rees (1995) quote a size of the 
order 3 x 10 8 m for a system that had an initial binary or- 
bital period of lOd and a companion of 15 solar masses. The 
IR limits for 4U0142+62 are similar and imply a somewhat 
smaller disk limit if the distance to the source is at the lower 
end of the proposed range of (i.e. 0.5 kpc) - maybe by up 
to a factor of 10 (i.e. a disk size of 7 x 10 r m). However, the 
upper bounds to the distance (1.5 kpc) give limits similar 
to the disk sizes of Podsiadlowski, Cannon and Rees (1995). 
Thus these new IR values are probing the expected disk sizes 
for a proposed TZO envelope. 



4 CONCLUSIONS 

The first deep infrared and optical images have been pre- 
sented of three proposed remnants of Thorne-Zytkow Ob- 
jects (TZO). In particular, the IR data go several infrared 
magnitudes deeper than previous observations and in at 
least one case, 1E2259+586, reveal the existence of weak 
objects within the X-ray error circle. It is argued, however, 
that in that case none of the objects seen is likely to be the 
binary companion to the X-ray source in that region. These 
data present new limits on any possible IR counterpart to 
the respective X-ray objects. 
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